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Abstract 

We report on the creation of new lasing modes with spatially nonuniform distributions of optical 
gain in a one-dimensional random structure. It is demonstrated numerically that even without 
gain saturation and mode competition, the spatial nonuniformity of gain can cause dramatic and 
complicated changes of lasing modes. New modes appear with frequencies in between those of the 
lasing modes with uniform gain. We examine some new lasing modes in detail and find they exhibit 
high output directionality. Our results show that the random lasing properties may be modified 
significantly without changing the underlying structures. 

PACS numbers: 42.55.Zz, 42.60.-v, 42.55.Ah 
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A conventional laser consists of a resonant cavity and amplifying material. The lasing 
modes have a nearly one-to-one correspondence with the resonant modes of the cold cavity 
Except for a slight frequency pulling, the lasing properties are usually determined by 
the cavities. Thus, cavity design is essential to obtain desirable lasing frequencies or output 
directionalities. The available lasing modes are typically fixed once the cavity is made. Finer 
control over lasing properties can be obtained, for example, by carefully placing the gain 
medium in a cavity to reduce the lasing threshold 2 ] or using specific pumping profiles to 
select lasing modes with desirable properties 0, 4, 5, 6]. However, once the laser cavity is 
made, it is very difficult to obtain new lasing modes that have no correspondence to the 
resonant modes of the cold cavity if nonlinearity is negligible. 

A random laser is made of disordered media and the lasing modes are determined by the 
random distribution of refractive index. Because of the randomness, it is difficult to inten- 
tionally produce lasing modes with desirable properties. To have more control over random 
aser properties, the structures themselves may be adjusted by selecting the scatterer size 
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11] and separation jl2|, |l3j] , changin g th e scattering structure with temperature 



15] or electric field [16], or creating defects [17] . For random lasers operating in the lo- 



calization region, spatially non-overlapping modes may be selected for lasing through local 
pumping of the random system [l{|. In the case of diffusive random lasers, far above the 
lasing threshold, nonlinear interaction between the light field and the gain medium alters 
the lasing modes [h]]. Without gain nonlinearity, local pumping and absorption in the un- 
pumped region can also change the lasing modes significantly [20] because the system size is 
effectively reduced. Recent experiments 21. 22] and numerical studies {^J show that even 
without absorption in the unpumped region, the spatial characteristics of lasing modes may 
vary with local pumping. In this case, the lasing modes still correspond to the resonant 
modes of the passive system. However, spatial inhomogeneity in the refractive index can 
introduce a linear coupling of resonant modes mediated by the polarization of gain medium 



24]. 



In this Letter, we demonstrate that new lasing modes can be created by nonuniform 
distributions of optical gain in one-dimensional (ID) random systems without absorption 
and nonlinearity. These new lasing modes do not correspond to the modes of the passive 
system or any lasing modes in the presence of uniform gain. They typically exist for specific 
gain distributions and disappear as the profile is further altered. They can lase independently 
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of other lasing modes when gain saturation is taken into account. The new lasing modes 
appear at various frequencies for many different gain distributions and can have highly 
directional output. These findings may offer an easy and fast way of dramatically changing 
the random laser properties without modifying the underlying structures. 

We consider a ID random system composed of N = 161 layers. Dielectric material with 
index of refraction ni = 1.05 separated by air gaps (712 = 1) resulting in a spatially modulated 
index of refraction n(x). The system is randomized by specifying different thicknesses for 
each of the layers as = (^1,2) (1 + vOi where (dx) = 0.1 /im and (c^) = 0.2 /zm are the 
average thicknesses of the layers, rj = 0.9 represents the degree of randomness, and ( is a 
random number in (-1,1). The length of the random structure L is normalized to (L) = 24.1 
/im. The index of refraction outside of the random media is n = 1. The above parameters 
give a localization length of £ ~ 240 fim at a vacuum wavelength A = 600 nm, which is the 
wavelength of interest in this work. 



The transfer matrix (TM) method developed in [23J is used to simulate lasing modes at 
the threshold with linear gain. A real wavenumber k = 2ir/\ describes the lasing frequency. 
Propagation of the electric field through the ID structure is calculated via the 2x2 transfer 
matrix M. The boundary conditions at the lasing threshold with only emission out of the 
system require M22 = 0. Linear gain is simulated by appending an imaginary part to the 
index of refraction n(x) = n(x) + mj/^x), where < 0. We neglect the change of the 
real part of the refractive index n(x) in the presence of gain. Spatial nonuniformity of gain 
is implemented by multiplying the imaginary part by a step function Je{x) = ©(— x + Iq), 
where x = is the left edge of the random structure and x = Iq is the location of the gain 
edge on the right side. 

Lasing frequencies and thresholds (gain required to induce lasing) are located by deter- 
mining which values of k and n^, respectively, satisfy M 2 2 = 0. Re[M 2 2] = (Im[M 2 2] = 0) 
forms real (imaginary) "zero lines" in the (k, n«) plane. The crossing of a real and imaginary 
zero line results in M22 = at that location, thus pinpointing a solution. We visualize these 
zero lines in Fig. [H by plotting — log 10 |Re M22I and — log 10 |ImM 2 2| together and using 
image processing techniques to enhance the contrast. We monitor the changes of such zero 
lines as the gain edge is moved gradually from Iq = L (uniform gain) within the wavelength 
range 500 nm < A < 750 nm. Appearances of new lasing modes with frequencies in between 
those of the lasing modes with uniform gain are observed. Figure [1] concentrates on a smaller 
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FIG. 1: (Color online) Left column: Real (green) and Imaginary (red) zero lines of M22 in the 
(k,rii) plane. Right column: Phase 9 of M22 in the (k,rii) plane. The horizontal axis is real 
wavenumber k in the range 10.48 /mi -1 and 10.64 /mi" 1 . The vertical axis is the imaginary part of 
refractive index rii in the range -0.0130 and -0.0326. The length Ig of the gain region is from top 
to bottom: (a-b) 14.961 /im, (c-d) 14.559 /mi, (f-g) 14.523 /mi, (h-i) 14.472 /mi, (j-k) 14.295 /mi. 

frequency range with Fig. Ufa) showing two of the lasing modes (marked mode 1 and mode 
2) resulting from the crossing of zero lines for Ig = 14.961 /im. Mode 1 at Ai = 592.6 nm has 
a lower lasing threshold than mode 2 at A2 = 598.0 nm. They are the only two lasing modes 
found within this frequency range and their frequencies are almost the same as those of the 
corresponding resonant modes of the passive system. As Ig is decreased, the lasing modes 
do not shift much in frequency. However, when Ig = 14.559 /im, the zero lines are joined as 
shown in Fig. [T](c). A new lasing mode, encircled in white, appears in between modes 1 and 
2. The spatial intensity distribution of the new lasing mode differs from those of modes 1 
and 2. As Ig decreases further, the zero lines forming mode 2 and the new mode pull apart. 
This causes the solutions to approach each other in the (fc, n^) plane, seen in Fig. [T](f), until 
becoming identical. The zero lines then separate resulting in the disappearance of mode 2 
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and the new mode as evidenced by Fig. QJh). The lines cross again for Iq = 14.295 /zm and 
the solutions reappear and move away from each other in the (k, rii) plane in Fig. [JQ). 

Verification of lasing mode solutions is provided by the phase of M 22 , calculated as 6 = 
atan2(Im[M 22 ], Re[M 22 ]). Locations of vanishing M 22 give rise to phase singularities. The 
phase change around a closed path surrounding a singularity is referred to as topological 



charge 25|, |26( . Two phase singularities are seen in Fig. (U^b) at the same locations as the 
zero line crossings in Fig. QJa). This verifies the authenticity of the lasing mode solutions. 
The phase singularity at the location of the new mode in the (k, n^) plane [Fig. QJd)] 
confirms that it is a genuine lasing mode in the presence of linear gain. The phase singularity 
associated with the new mode is of opposite charge to the existing ones. As Iq is reduced, 
two oppositely charged phase singularities move closer [Fig. QJg)] and eventually annihilate 
each other at Iq = 14.295 /im [Fig. QJi)]. As already mentioned, this process reverses itself 
and the two lasing modes reappear in Fig. []Jk). 

For a more thorough study of the new lasing modes and confirmation of their existence 
in the presence of gain saturation, we switch to a more realistic gain model including non- 
linearity. The Bloch equations for the density of states of two-level atoms [27(] are solved 
together with the Maxwell's equations with the finite-difference time-domain method [281 ] . 
The phenomenological decay times due to the excited state's lifetime T\ and decoherence T 2 
are included. The gain spectral width is given by Au; a = 1/Ti + 2/T 2 We also include 
incoherent pumping of atoms. The rate of atoms being pumped from the ground state to 
the excited state is proportional to the ground state population, and the proportional coeffi- 
cient P r is called the pumping rate. The resulting Maxwell-Bloch (MB) equations are solved 
numerically with the spatial grid step Ax = 1.0 nm and the temporal step At = 3.3 x 10~ 18 
s. The atomic density N at0 m/V = 1.8 x 10 13 cm -3 . Nonuniform gain is simulated by having 
the two- level atoms only in the region < x < Iq. 

By setting the atomic transition wavelength A a to coincide with the wavelength of mode 
1, 2 or the new mode and using a narrow gain spectrum, we are able to investigate the three 
lasing modes separately. Au a is chosen to be less than the mode spacing to ensure single 
mode lasing (at smaller pumping rates). At Iq = 14.295 /im, the wavelength difference 
between mode 2 and the new mode, which is smaller than that between mode 1 and the 
new mode, is A 2 - \ nm = 1.4 nm. We set T x = 1.0 x 10~ 12 s and T 2 = 0.73 x 10" 12 s so that 
the gain spectral width in terms of wavelength is AA a = 0.71 nm. Initially all atoms are 
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FIG. 2: (Color online) Steady-state output intensity vs. pumping rate P T from Maxwell-Bloch 
simulations with different gain spectra. The wavelength of the atomic transition A a is equal to Ax 
(red crosses), A2 (blue open diamonds), and X nm (black circles). 

in their ground state and the system is excited by a Gaussian-sinusoidal pulse with center 
wavelength A = A a and spectral width AA = AA a . When the pumping rate P r is above a 
threshold value, the electromagnetic fields build up inside the system until a steady state is 
reached. 

Figure [2] shows the steady-state output intensity with A a = Ai, A 2 , or A nm as P r is varied. 
P r = 1/Ti corresponds to the transparency point, namely, the excited state population of 
atoms is equal to that of the ground state. The lasing threshold pumping rate for mode 1 is 
reached first at P r = 1.9/Xi, then mode 2 at P r = 2.1/Xi and the new mode at P r = 3.0/Ti. 
These thresholds agree qualitatively with the values of the TM calculation. 

When A a = A nm , the first lasing mode is the new mode, instead of mode 1 or 2. Figure 
El^a) shows the output emission spectrum just above the lasing threshold at P r = 3.0/Ti. It 
consists of a single lasing mode with the wavelength equal to that of the new mode calculated 
with the TM method. The spatial intensity distribution obtained from the Maxwell-Bloch 
(MB) calculation is compared to that from the TM calculation in Fig. [3]^b). The MB 
distribution \4>mb(x)\ 2 is found by integrating the intensity over one optical period. It is 
then normalized to the TM distribution \4>tm(x)\ 2 as Jq \4>MB(x)\ 2 dx = \(f)TM{x)\ 2 dx. 
The two intensity distributions are almost identical. The average percent difference between 
them is 7.77%. This result indicates the nonlinear effect due to gain saturation is small 
when the pumping rate is just above the lasing threshold. When the peak of the gain 
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FIG. 3: (Color online) (a) Emission spectra taken from two Maxwell-Bloch simulations with A a 
equal to X nm (solid line) and A2 (dashed line). In both cases, AA a = 0.71 nm, Iq = 14.295 //m. 
P r = 3.0/Ti for the solid curve and 2.1/Ti for the dashed curve. The emission intensity reaches a 
steady state by a simulation time of t = 300 x 10~ 12 s. The spectra are taken from 1400 x 10~ 12 
s to 1700 x 10~ 12 s. The MB intensity distribution \4>mb{x)\ 2 (black dashed line) is compared to 
the TM intensity distribution \(Ptm(x)\ 2 (red solid line) for the new lasing mode (b) and mode 2 
(c). The inset in (b) is an expansion of the curves for fim < x < 10 fim. 

spectrum is shifted from A nm to Ai or A 2 , the first lasing mode is switched to mode 1 or 2. 
Figure [3](c) plots the spatial intensity distribution of mode 2 obtained by the MB calculation 
with A a = A2 and P r = 2.1/Ti as well as that obtained by the TM calculation. The two 
distributions are almost the same and they are different from the distribution with uniform 
gain. Comparing Fig. Wh) to (c), we see the spatial intensity distribution of the new lasing 
mode differs significantly from that of mode 2 within the gain region. Outside the gain 
region the two distributions are not much different because their wavelengths are very close. 

When optical gain is located on the left side of the structure, we observe that the intensity 
distributions of new lasing modes are heavily concentrated on the right side of the gain 
region. This makes the emission intensity through the right boundary of the random system 
much larger than that through the left boundary. We calculate the the ratio of right to left 



7 



output flux S = \4>{x = L)\ 2 /\(f>(x = 0)| 2 . For the new lasing mode in Fig. Mjo) S nm = 40, 
indicating the laser output is mostly to the right. As a comparison, S\ = 1.1 for mode 1 and 
S 2 = 3.3 for mode 2. Thus, the new lasing mode has much more directional output than 
modes 1 and 2. 

Because of the excellent agreement found between the MB and TM calculations, we 
conclude that new lasing modes do appear in random lasers with spatially nonuniform 
distributions of optical gain. Typically, as in the case studied here, they are sensitive to 
the spatial gain distribution and disappear if the distribution is altered slightly. These 
new lasing modes offer more control of random laser performance as their properties such 
as frequency and output directionality can be quite different from those of existing lasing 
modes. Moreover, the properties of new lasing modes can be easily altered by varying the 
spatial profile of the pump beam, without modifying the random structures. 
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